Dual neutral variables and knot solitons in triplet superconductors 
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In this paper we derive a dual presentation of free energy functional for spin-triplet superconduc- 
tors in terms of gauge invariant variables. The resulting equivalent model in ferromagnetic phase has 
a form of a version of the Faddeev model. This allows one in particular to conclude that spin-triplet 
superconductors allow formation of stable finite- length closed vortices (the knotted solitons). 



In this paper we discuss spin-triplet superconductors 
and show that under certain conditions these supercon- 
ductors should allow formation of stable finite-length soli- 
tons. The discussions of the stable finite-length solitons 
in 3-f 1 dimensions have a long history This con- 

cept was first discussed in mathematical physics: In 
L.D. Faddeev introduced a version of nonlinear 0(3) a- 
model which allows formation of the stable finite-length 
solitons which may have a form of a knot or a vortex 
loop. The stability of these defects in the Faddeev model 
is ensured by a higher-order derivative term: 



Ht. 



dr 



a|5fcn|2 + — (n-a.nx 9,n)2 



(1) 



where n is a three-component unit vector. 

Recently it was realized that versions of this model 
are potentially relevant for description of many differ- 
ent physical systems ranging from infrared limit of QCD 
to superconductivity in transition metals Rele- 
vance of this model for condensed matter physics was 
pointed out in Ref. were it was found that in two- 
band superconductors there exists a mapping between 
two-flavour Ginzburg-Landau (GL) functional and a ver- 
sion of 0(3)-symmetric Faddeev model. The knotted 
solitons are much more complex and structurally com- 
plicated topological defects than Abrikosov vortices and 
thus its realization in superconductors should open ex- 
ceptionally wide range of possibilities of studies of var- 
ious phenomena associated with them. The remarkable 
circumstance is that the studies of the properties of these 
defects in superconductors may result in a "feedback" for 
the discussions of properties of possible similar defects in 
the infrared limit of QCD where it was argued that it 
describes glueballs 

So far the two-band superconductors were the only 
known condensed matter system which is described by 
a version of the Faddeev model and allows formation of 
the knotted solitons In this paper we discuss pre- 
sentation of the Ginzburg-Landau functional for a single- 
condensate spin-triplet superconductor in gauge- invariant 
variables by means of an exact duality mapping which 



also leads to another version of the Faddeev model in 
ferromagnetic state of the condensate. A remarkable 
fact which we observe below is that the resulting effec- 
tive model for ferromagnetic spin-triplet superconductor 
is surprisingly similar to the model discussed in albeit 
the physical origin of this model in spin triplet supercon- 
ductors is principally different from the model in . 

We would like to emphasis that in a charged spin-1 
Bose condensate the stable knotted solitons (which en- 
ergy is a nonmonotonic function of its size) are the coun- 
terpart of the Volovik-Mineev vortices characterized by a 
nontrivial Hopf invariant in a neutral spin-1 superfluid 
(discussed in the pioneering papers on superfluid '^He 
|, |l^, 0). The defects considered in |, |l^, |ri|, |l| 
have similar topology but its energy is proportional to its 
size and thus these defects are unstable against shrinkage. 

We begin with reminding a reader some important fea- 
tures of the neutral spin-1 condensate following the paper 
jlSf . This will be particularly helpful in the discussion 
of effects of coupling to gauge field and differences of the 
knotted solitons and Volovik-Mineev vortices. 

We write the order parameter of the spin-1 Bose con- 
densate as 5'a(r) = ^/n{x)(ai^) where (a = 1, 0, —1) with 
C being a normalized spinor • ^ = 1. Then the energy 
functional for a neutral spin-1 system can be written as 
ill: 



K = 



[ dr— (V^A^)2 + —n(VCY 
J 2M^ ^ ' 2M ^ ^' 



y [co -f C2 < F >2] 



/in 



(2) 



where < F Ca^abCb is spin. Degenerate spinors 
are related to each other by gauge transformation e'^ 



and spin rotations lA (a, /3, r)=e 



where 
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(a, (3, t) are the Euler angles. The ground state structure 
of (r) can be found by minimizing the energy with 
fixed particle number p^ . Below we shall be mainly in- 
terested in ferromagnetic state. This state emerges when 
C2 < 0. The energy is minimized by < F 1 and the 
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ground state spinor and density are ul 







• 









e-'"cos2| 
/2cosfsinf 
e^^sin^f 



n°(r) = 



1 



Co + C2 



(3) 



From this equation it is seen that in ferromagnetic case 
there exists an equivalence between gauge transformation 
9 and spin rotations r so the symmetry group is 50(3) 
^ . Given the expression for the ground state spinor (||) 
one can immediately derive the superfluid velocity for a 
neutral spin-1 ferromagnetic Bose system |l3|: 



[compare with two-gap superconductor |^]: 

c Tic 
Bfc = - ^ [V,C, ~ V,C,] + ^ (s • V,s X V,s), (7) 

where = ^ and we introduced the following nota- 
tions: 



8= (sin /3 cos a, sin /3 sin a, cos /3); 



en 



(8) 



(9) 



Let us rearrange terms in the Ginzburg-Landau func- 
tional (H). First let us rewrite the second term in (||) as 
follows: 



M 



[Vie 



cos /3Va] 



(4) 



Let us now turn to a charged spin-1 Bose condensate. 
We have the following expression for the free energy of 
spin-1 superconductor in the "ferromagnetic" state: 



F = dr 



fc2 j-2 

n n n 
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(5) 



with ground state spinor and density being given by (|^). 
Consequently the equation for the supercurrent is: 



ihen 
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Mc 
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r) - cos/3Va r^— A 

Mc 



(6) 



From this equation it is clearly seen that the supercur- 
rent depends not only on phase gradients but also on spin 
texture what is in direct analogy with the situation in 
neutral case (^). However, the properties of the charged 
spin-1 condensate (in particular topological defects al- 
lowed by the system) are principally different from the 
neutral case that can be seen if we eliminate gauge field 
by a duality transformation to neutral variables which ex- 
plicitly shows physical degrees of freedom in the system. 
Such a mapping, which is actually different comparing to 
the problem in Ref . Q| , is discussed below. 

We should emphasis that in the charged case (||) the 
free energy features a contribution from magnetic field 
B^/Stt which can be external or self- induced or both. 
We also stress that we consider below a superconductor 
in a simply-connected space, that is, our defects do not 
feature zeroes of the order parameter. In this case taking 
curl from both sides of (^), (and taking into account that 
in a simply-connected space for a regular function holds 
identity curl\/{9 — t) =0) we arrive at the following 
equation for the magnetic field in triplet superconductor 
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where 



j = (zCVCa - *CaVC) . 

Then we observe the following circumstances: 



(10) 



(11) 



-(V(0- 



COS (3 



COS/3 



(V(0-r)) 



2 sin P 



From the above expression it follows that: 



(Va)2 = -(Vs)2(12) 



Now we can express the Ginzburg-Landau functional for 
spin-1 charged ferromagnetic Bose condensate in the form 
involving only gauge-invariant variables (compare with 
[H ). That is, with this procedure we ehminate the gauge 
field A and the variable — r in favor of the gauge- 
invariant variables s and C. With it the free energy func- 
tional for the ferromagnetic p-wave superconductor be- 
comes: 



F 



dr 



2M 



(vVH)2 
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(13) 
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The above expression is a version of the Faddeev nonhn- 
ear 0(3) cr-model (|l|) introduced in The remarkable 
circumstance is that this effective model for a spin-triplet 
superconductor is very similar to the effective model for 
two-band superconductors where Cooper pairs have spin- 
but the system possesses a hidden 0(3)-symmetry |^ 
and can be described by means of a three-component 
unit "pseudospin" vector. The resulting model explicitly 
displays only physically relevant degrees of freedom and 
indeed does not depend on 0, r and A. The new vari- 
ables are: n = I'I'P, the vector s which position on the 
unit sphere is characterized by the angles a and (3 
and the massive vector field C. The effective action 
contains fourth-order derivative term for the vector field 
s. The field s is also coupled to the massive vector field 
C. Thus we can see that the properties of the charged 
spin-1 condensate are very different form the neutral con- 
densate (H). Here we would like to stress that the proce- 
dure with which we arrived to (|l^) is the exact change of 
variables with no approximations were used which means 
that there exists the exact mapping between the models 
( p^ ) and (^). The effective model (|l^) displays features 
of spin-triplet superconductors which can be easily over- 
looked in the presentation of the effective functional in 
the form (||). 

In particular, from (^3|) we can conclude that the sys- 
tem allows stable knotted solitons characterized by a non- 
trivial Hopf invariant which stability is ensured by the 
fourth-order derivative term (the Faddeev term). This 
term is not explicitly present in the free energy func- 
tional in the Ginzburg-Landau form (j^). It stems from 
the term B^/Stt and has the physical meaning of the self- 
induced magnetic field in a presence of a nontrivial spin 
texture (in this paper we discuss a system with no ap- 
plied external fields) . Indeed this effect is possible due to 
the feature of the ferromagnetic state of p-wave conden- 
sate where the magnetic field may be induced by a spin 
texture. 

Since the effective models for spin-triplet superconduc- 
tor and two-band superconductor appear being formally 
very similar we refer a reader to the paper for a de- 
tailed description of the knotted solitons in the model 
( p^ ) , whereas below we outline differences of the proper- 
ties of solitons in these systems. One of the differences is 
that in the spin-triplet case the knot soliton may form as 
a nontrivial texture with no inhomogeneities in Cooper 
pair density. In particular that means that there is no 
mass of the components of the vector s whereas in two- 
band case [|j the 123 component of the unit vector n is re- 
lated to relative local Cooper pair densities and is indeed 
massive. In principle, in triplet case, the components of 
the field s may also acquire a mass if to take into account 
spin-orbit interaction which would give s an energetically 
preferred direction, here we assume spin-orbit interaction 
to be small comparing with mass for the field C. Under 
the assumption of small spin-orbit interaction the char- 



acteristic size of a knotted soliton is determined by a 
competition of the second-order and fourth-order deriva- 
tive terms for s and thus it is of order of magnitude of 
magnetic field penetration length JlSt . 

We stress that in two- band superconductor in the 
points where the unit vector n is situated on the south or 
north poles of unit sphere S"^, the densities of Cooper pair 
of flavour 1 or flavour 2 vanish correspondingly. Since a 
knot soliton is characterized by a Hopf invariant that 
means that the vector n necessary hits poles of the unit 
sphere which necessary results in zeroes of Cooper pairs 
densities. Also the contribution to the self-induced mag- 
netic field associated with the term n ■ diH x djii van- 
ishes in these areas. In contrast, the knotted soliton 
with the same topology in spin-triplet superconductor is 
a nontrivial configuration of spin texture and it does not 
feature zeroes of density of Cooper pairs. When spin- 
orbit interaction is small we can choose e.g. that the 
vector s assumes at infinity the value corresponding to 
the north pole of the unit sphere and in the knot soli- 
ton core reaches the south pole. The configuration of the 
self-induced magnetic field associated with contribution 
of s • diS X djS in spin-triplet case is following: the mag- 
netic field vanishes indeed at spatial infinity (north pole 
of the unit sphere) and also the magnetic field vanishes 
in the knotted soliton "core" In between the core 

and "knot boundary" the self-induced magnetic field as- 
sociated with the contribution of s • diSX djS has a helical 
geometry characterized by a corresponding Hopf invari- 
ant. 

Thus the knot soliton in triplet superconductor has dif- 
ferent structural features and different underlying physics 
than the knotted solitons in two-gap superconductor al- 
beit the main feature that these defects have in common 
is the self-induced nontrivial configuration of magnetic 
field which amounts to a Faddeev term in the effective 
action and which stabilizes the size of the soliton. In 
other words a shrinkage of a knotted soliton would re- 
sult in increasing energy of self-induced magnetic field 
which gives this defect energetic stability [compare with 

011]. 

Let us also emphasis that the "polar" phase of triplet 
superconductors is not described by a version of the Fad- 
deev model and does not allow formation of knotted soli- 
tons in contrast to the ferromagnetic phase. The prop- 
erties of a neutral polar phase were investigated in [ p^ . 
This state appears in the case when C2 > 0. The energy 
is minimized by < F 0. The spinor ( in the ground 
state was calculated in ||l^. In a charged counterpart of 
the system in the polar phase considered in |l3|| , it can 
be easily seen that the supercurrcnt does not depend on 
the spin texture and this formally results in the absence 
of a Faddeev term in the effective action. 

In conclusion, we derived an exact equivalent presenta- 
tion of the free energy functional for ferromagnetic spin- 
triplet superconductor. The derived functional in dual 
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gauge-invariant variables explicitly displays the physical 
degrees of freedom in the system. In particular it al- 
lows us to conclude that the ferromagnetic spin-triplet 
superconductors allow formation of stable knotted soli- 
tons. This is in contrast to a neutral spin-1 condensate 
where the finite-length defects characterized by a non- 
trivial Hopf invariant arc not energetically stable [p[-p2[. 



The amazing fact is that the derived model is very simi- 
lar to the model considered in |2| despite it was derived 
from a very different Ginzburg-Landau functional. It 
shows that this version of the Faddeev model is a rather 
generic model for various superconductors. The interest- 
ing question is an experimental procedure how to cre- 
ate and observe the knotted solitons. Indeed in contrast 
to Abrikosov vortices these defects can not be created 
by simply applying an external magnetic field. On the 
other hand the situation is simplified by the key feature 
of these defects, which is, once they are formed they are 
stable being protected against decay by an energy bar- 
rier. One may expect that e.g. a rapid cooling of a 
system from above to below critical temperature in e.g. 
applied random fields should indeed result in a formation 
of a certian density of topological defects. Apparently an 
ensemble of these defects, because of its complex struc- 
ture, and very different nature comparing to Abrikosov 
vortices, should exhibit many unconventional phenom- 
ena and thus should be a very interesting object for ex- 
perimental studies, especially since the triplet supercon- 
ductivity has been established experimenally in several 
compounds (e.g. UPt^ and Sr2RuOi jl^). Finally, we 
would like to especially stress that it was suggested in 
that knotted solitons could play an important role in 
the infrared limit of QCD. The fact that similar defects 
should be present in p-wave and two-band superconduc- 
tors and the macroscopic quantum origin of the topologi- 
cal defects in condensed matter systems which implies its 
rather direct observability means that the spin-triplet su- 
perconductors along with two-band superconductors may 
in some sense serve as "a testing laboratory" for the in- 
frared limit of QCD li- 
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The coupling to the field C also should affect a knot 
soliton: a nontrivial configuration of texture s induces 
a magnetic field which in turn should induce screening 
Meissner current which amounts to a nontrivial configu- 
ration of C. However the coupling to C could not kill the 
knot soliton. Indeed let us assume a knot is shrinking to a 
characteristic size smaller than the magnetic field pene- 
tration length A, then we have a nontrivial configuration 
of s which amounts to self-induced magnetic field con- 
figuration proportional to s • ViS x VjS which becomes 
more and more singular if the solition is shrinking. The 
direction of this self-induced field also changes inside the 
knot since the configuration of s is helical and is char- 
acterized by a Hopf invariant. On the other hand the 
Meissner current C screens magnetic field in a supercon- 
ductor over a characteristic length scale A thus if a tex- 
ture s shrinks at characteristic sizes smaller than A we 
have a fast variable s and a slow variable C over which 
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we can average and thus at the length scales smaller than 
A the model reduced to the model with an arbitrary 
accuracy thus the knot solitons in the model (^3|) are 
stable. This can be also understood intuitively: at the 
length scales smaller than A the system almost does not 
display Meissner effect and thus the screening current 
cannot affect the highly nontrivial configuration of self- 



induced magnetic field proportional to s ■ V^s x VjS. 
[19] We should emphasis that under the core of the knotted 
soliton we understand here the closed line (which may 
form a loop or a knot) where the vector s reaches the 
south pole. There are no zeros of Cooper pair density in 
what we call here the "core". 



